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Abstract 

We investigate the gradual emergence of the disorder-related 
phenomena in intermediate regimes between a deterministic 
periodic Bragg grating and a fully random grating and high¬ 
light two critical properties of partially disordered Bragg grat¬ 
ings. First, the integral of the logarithm of the transmittance 
over the reciprocal wavevector space is a conserved quantity. 
Therefore, adding disorder merely redistributes the transmit¬ 
tance over the reciprocal space. Second, for any amount of 
disorder, the average transmittance decays exponentially with 
the number of grating layers in the simple form of exp (—rjN) 
for sufficiently large N, where 77 is a constant and N is the 
number of layers. Conversely, the simple exponential decay 
form does not hold for small N except for a highly disordered 
system. Implications of these findings are demonstrated. 

I. Introduction 

An optical Bragg grating is created by periodic variations in 
the permittivity of dielectric material. Depending on the size 
of the periodicity, the grating can act as a highly reflective 
(bandgap) or a highly transmittive (bandpass) mirror (lj. A 
bandgap or bandpass is not limited to an optical Bragg grating: 
in solid-state physics, a periodic lattice of atoms or molecules 
in a crystal modulates the background potential for electron 
wave function and similar badgap and bandpass characteristics 
are observed. Understanding the electronic band structure is 
at the foundation of solid-state physics. 

In practice, no potential is perfectly periodic and no re¬ 
fractive index modulation is disorder-free. Therefore, it is 
important to investigate the impact of disorder on the band 
structure in electronic systems, optical gratings, and other 
coherent wave systems [ 2 ). There have been numerous studies 
of disordered coherent wave systems, much of which follows 
the seminal work of Philip Anderson ( 3 ]]. He showed that for a 
sufficiently high level of disorder the electronic wavefunction 
is exponentially attenuated, resulting in diminished transport 
and conductivity 0, as also confirmed in subsequent work by 
others | 5 ]|—fl6). 

A one-dimensional (ID) periodic Bragg grating is simple 
and elegant-its band structure is easy to understand and 
analytical calculations can often be performed to study its 
characteristics GZ) Moreover, many of its properties are 
generic and can be extended to 2 D and 3 D gratings. A 
highly disordered ID grating made from a random stack of 
dielectric slabs has also been treated analytically by Berry and 


Klein p8| . They showed that the average transmitted intensity 
drops exponentially with the number of slabs. 

The problem of a perfectly periodic Bragg grating is well 
understood, and so is that of a highly disordered grating. The 
intention of this article is to shed further light on intermediate 
situations with moderate disorder fT^-pO]]- We explore the 
transition from a deterministic periodic Bragg grating to a 
fully random grating, and highlight the gradual emergence 
of the disorder-related phenomena in intermediate regimes. 
In particular, we present two critical properties of partially 
disordered ID Bragg gratings. Because our analysis is carried 
out for a stack of dielectric slabs similar to Ref. [ 18 ], we 
present these properties in the language of dielectric slabs for 
simplicity: 

First we find a globally conserved quantity in the reciprocal 
wavevector space that is quite instrumental in visualizing 
the wave localization behavior in a partially disordered stack 
of dielectric slabs. We show that the total area under the 
curve of the logarithm of the transmittance plotted in the 
reciprocal wavevector space is the same for all Bragg gratings, 
regardless of the amount of disorder. As a consequence of this 
conservation law, if the transmittance is higher in some region 
of the reciprocal space, it has to be lower in another region to 
conserve the total area. For example, the presence of a strong 
bandgap for a periodic Bragg grating necessitates a strong 
bandpass in a different region of the reciprocal space. 

Second we observe that for any amount of disorder, the 
average transmitted intensity decays exponentially with the 
number of slabs as exp(— rjN) for sufficiently large N. r\ is a 
constant whose value depends on the amount of disorder, and 
N is the number of slabs. Even in the absence of disorder, the 
simple exponential decay law holds in the bandgap region for 
large N. Conversely, the simple exponential decay form does 
not hold for small N except for a highly disordered system. 

In the following, these ideas will be presented in the frame¬ 
work of optical transmission through a ID stack of dielectric 
slabs. The slabs are assumed to be identical with refractive 
index ri2 and are embedded in a background dielectric of 
refractive index n\. The disorder is introduced by randomizing 
the location of the slabs. 

II. Background 

The transmission (transfer) matrix for a lossless reciprocal 
mirror (dielectric slab) can be expressed as 
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t and r are the amplitude transmittance and reflectance, 
respectively, and are functions of the dielectric constants and 
geometric properties of the slab and frequency and incidence 
angle of light JTJ . They also satisfy the losslessness condition 
|t| 2 + |r| 2 = 1. The latter simplified form is obtained by 
formally incorporating a thin layer of the background rt\ 
dielectric in the n 2 slab, where r = \t\ and p = |r|. The 
latter form is used in the following discussions as it does not 
affect the generalities of the presented arguments. 


The transfer matrix of an array of N identical dielectric 
slabs with varying separations (gaps) of the background di¬ 
electric material can be expressed as 
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where J is the translation matrix accounting for the phase 
accumulation (p n in the nth gap. Here, is assumed to 

be the identity matrix. The total wave transfer matrix of the 
system is given by 
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where the total intensity transmittance is given by the (2,2) 
element of the T(N,t; {</?„}) = |M^ ) |“ 2 - 


III. Transmission through a periodic vs. random 
stack 



Fig. 1. Transmittance is plotted as a function of the gap phase for a 20-layer 
periodic Bragg grating with r = 0.8. 



Fig. 2. Average transmittance is plotted in logarithmic scale as a function 
of the average gap phase. Each grating has N = 50 layers with r = 0.5. 
Different curves are for different values of the disorder parameter a. The area 
under each curves is the same. 


The most common application of the preceding analysis is 
to study periodic Bragg gratings. For a periodic grating where 
the gap thickness between consecutive slabs is identical, the 
accumulated phase in each gap takes a common value (p. In 
Fig. [T] the total transmittance T is plotted for a periodic stack 
of N = 20 dielectric mirrors with r = 0.8 as a function (p, 
where one can clearly see the familiar bandpass and bandgap 
regions. 

In practice, it is impossible to maintain a uniform gap (fixed 
(p) between all slabs, and some randomness is inevitable. A 
convenient way to parametrize a partially disordered grating is 
to assume that the accumulated phase in each gap is a random 
number with an average value of (p with some probability 
for variation around the average. We adopt the definition in 
Eq. |4j which states that (p n is chosen from a uniform random 
distribution between <p — an and (p + an. The disorder level 
is parametrized by a: the deterministic periodic grating is 
identified by a = 0 and a fully random grating by a = 1. 

= ( P + ^n) 0 n (E a x unif[—7r, 7r], 0 < a < 1. (4) 

The transmittance for different levels of randomness char¬ 
acterized by a is plotted in Fig. [2] as a function of (p. Here, 
N = 50 and r = 0.5 is assumed and the vertical scale 
is logarithmic. For each value of a , the transmittance curve 
is properly averaged (proper averaging will be explained 
shortly), for an ensemble of 10,000 different gratings, each 
for a different set of {p n } values. 

For a periodic Bragg grating with a = 0, an analytical 


formula exists for transmittance: 
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The transmittance curve in Fig. [l] and the black solid curve 
in Fig. [2] follow Eq. [5] The oscillatory bandpass regions 
correspond to | cos ((p) \ < r where is real, while the bandgap 
region correspond to | cos(</?)| > r, where <f> is imaginary. In 
the bandgap region, it is more convenient to redefine n and 
T> as 


T sinh(7VT>) , . /cos(<p)\ 

*N= . $ = cosh 1 -^ . (6) 

smh(<3>) V r J 

For a partially random grating with a > 0, the transmittance 
curve would look rather noisy. The smooth transmittance 
curves in Fig. [2] are only obtained after properly averaging over 
many gratings. Because the transmittance is a multiplicative 
quantity, proper averaging of the random transmitted intensity 
is done by averaging the logarithm of the transmittance 0. 
fl8] l. Using Eq. [3j we can formally express the averaging as 

N n(p-\-Ot 7T 7 

iog ( T "“) = (n J f _„ 2^) ioe ( r «» , "») 

= 27V log(r) + U(7V, r, (p, a), (7) 


where 2N log(r) comes from the overall multiplicative factor 
in Eq. [3] and fi(r, N,{(p n }) comes from averaging the matrix 
multiplication part. 




























































For the case of a totally random dielectric stack with 
a = 1, Berry and Klein p~ 8 | have rigorously shown that 
Q(N,r,(p,l) = 0 ; therefore, T avg = exp(27V log(r)) and is 
independent of the value of (p. The numerical simulation of 
Fig. [2] in red dotted line agrees with the analytical results, 
where log(7I V g) = 2A/Tog(r) = —69.31 for N = 50 and 
r = 0 . 5 . 

For intermediate values of a between 0 and 1, a closed- 
form solution is not available for T avg and one must resort to 
numerical plots similar to those presented in Fig. [2] However, 
all these curves follow a remarkable conservation law, which 
is a consequence of dp Q = 0. It can be verified that by 
changing the amount of randomness via a, the averaged trans¬ 
mittance profile redistributes itself over the mean accumulated 
gap phase p in a such way that it preserves the area under the 
logarithm-transmittance curves in Fig. [2j 

■A = J dp log (r avg (iV, r, p, a)) = 27rA r log(r). (8) 

The area is independent of the randomness level a and is 
the same for all curves. This conserved transmittance area 
theorem gives special status to the case of a totally random 
dielectric stack marked with a = 1 and the conventional 
periodic Bragg grating marked with a = 0. The totally 
disordered grating incorporates all possible values of phase 
and p loses its special standing; therefore, the transmittance 
becomes independent of p by uniformly spreading the avail¬ 
able conserved area A over all values of p. Conversely, 
the periodic Bragg grating with a = 0 provides the most 
nonuniform distribution of the available conserved area A over 
the space of p with highly depressed values of transmittance 
in the bandgap, accompanied by a large transmittance in the 
bandpass to compensate. 

In Appendix, we offer a proof of the conserved area 
theorem. It is shown that the conserved area theorem applies 
to the transmittance for a general partially random grating 
and averaging over {6 n } is not required ({6 n } was defined in 
Eq- 0 . In the absence of ensemble averaging, the transmittance 
curve would look rather noisy but still satisfies the conserved 
area theorem. 

The conserved area theorem provides a very intuitive and 
useful approach to visualize the impact of partial disorder. It 
is important to note that the accumulated phase values in the 
gaps are proportional to the wavevector k , where p n = n\kd n 
and d n is the random thickness of the nth gap. Each curve in 
Fig. [2] should be regarded as corresponding to a grating with 
an average gap thickness J, where p = n\kd. Therefore, the 
horizontal axis p in Fig. [2] actually represents the wavevector, 
and the conserved area theorem is a statement about the 
integral of the average of the logarithm of the transmittance 
over the reciprocal space. 

IV. Transmittance scaling with the number of 
mirrors 

The exponential decay of transmittance (relative optical 
intensity) is regarded as the main signature of Anderson 
localization. For the case of a totally random stack marked 
with a = 1 , it was previously shown that the average trans¬ 
mittance scales like T avg = exp (—77 AT), where the exponent 


(p = 0 



Fig. 3. Average scaled decay exponent 77 is plotted as a function of the 
number of mirrors N for different values of the disorder parameter ol.t — 0.5 
and (p = 0 are assumed here. 77 is meaningful only at integer values of N, 
but the points are connected for better visual clarity. 

is proportional to the number of mirrors. The scaled exponent 
r] is independent of AT and is given by r] = —21og(r) > 0. 

In this section, it is argued that neither the exponential decay 
nor its universal scaling with AT is specific to the case of 
a totally random stack marked with a = 1. Rather, for a 
sufficiently large number of mirrors AT, T avg always scales 
like exp(— rjN) for all gratings. The only exception is for the 
periodic Bragg grating marked with a = 0 in the bandpass 
region where | cos(<^)| < r. 

We start by studying the periodic Bragg grating with a = 0. 
According to Eq.[5j the transmittance in the bandpass is a non¬ 
decaying oscillatory function of AT. However, in the bandgap, 
it is an exponentially decaying function of AT for large AT. This 
can been seen by using Eq. [ 6 ] and noting that sinh(AT<F) ^ 
exp(AT<f>)/2 for large AT. After a few simple algebraic steps, 
it can be shown that log (T) ~ 2AT<f>; therefore, 77 = — 2<f> for 
large AT. For AT = 1, Eq. [5] shows that log(T) = 21og(r). In 
Fig- [3] r] is plotted as a function of AT for different values of 
a, all for p = 0 and r = 0.5. The black solid line corresponds 
to a = 0, where 77 starts at 21og(r) « —1.37 for AT = 1 and 
saturates at 77 = 2cosh _ 1 (l/r) « 2.63 for large AT. 

The plots of 77 as a function of AT for different values of 
a follow a similar behavior to that of the periodic grating 
explained above. They all start at 77 = 21og(r) for N = 1 
(see Eq. [5]) and asymptotically approach a limit for large N. 
Therefore, the average transmittance decays exponentially for 
all cases at large N. In Fig. [3j which is specific to (p = 0, the 
largest asymptotic value for 77 is obtained for a = 0. This is 
not surprising, because (p — 0 corresponds to the bottom of 
the bandgap of the periodic Bragg grating in Fig. [2] 

A similar behavior is observed in Fig. |4j which is the same 
as Fig. [3j except (p = 7 r/ 2 . The main difference is for the 
solid black line, which relates to the periodic Bragg grating 
with a = 0 : 77 oscillates and asymptotically approaches zero 
as 1/AT; therefore, transmittance does not decay exponentially 
even for large N. This behavior is expected because (p — n/2 
is at the center of the bandpass for the Bragg grating. When a 
slight amount of disorder marked by a = 0.1 is introduced, the 
77 -curve (long-dashed blue) goes through a brief oscillation and 
then approaches asymptotically to a nonzero value, indicating 
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Fig. 4. Similar to Fig. [5] except ip = 7 t/2 . 


that for large N the 7^ V g = exp (— 777 V) behavior is observed. 
Increasing randomness to a = 0.5 increases the scaled decay 
exponent monotonically. Increasing a beyond 0.5 initially 
increases 77 , but eventually rolls it back to make the 77 -curve 
for a = 1.0 the same as that of a = 0.5. 

The universal exponential decay behavior at large N is 
expected because all ID disordered systems are Anderson 
localized [4]. However, the converse observation that for small 
N , the transmittance curve can significantly deviate from 
the simple exponential decay form despite being Anderson- 
localized is worthy of special attention. The exponential decay 
is regarded as the main signature for Anderson localization 
and the stated converse observation signifies the importance 
of asymptotic analysis in establishing Anderson localization. 
In other words, large deviations from an exponential decay in 
the first few layers do not exclude the possibility of Anderson 
localization. 


V. Conclusions 

The transition from a deterministic periodic Bragg grating 
to a fully random grating and the gradual emergence of the 
disorder-related phenomena in the intermediate regimes is 
studied in detail. It is shown that for any amount of disorder the 
average transmittance decays exponentially with the number of 
grating layers in the simple form of exp (—rjN) for sufficiently 
large N, where 77 is a constant. For a highly disordered grating, 
this simple form is true regardless of the number of mirrors N. 
Conversely, for small N, the transmittance curve can deviate 
from the simple exponential decay form and can even oscillate 
as it decays. Whether increasing the disorder increases or 
decreases the decay exponent depends on the location of the 
grating in the reciprocal wavevector space, the amount of 
disorder, and other parameters defining the grating. 

It is also shown that the integral of the logarithm of 
the transmittance over the reciprocal wavevector space is a 
conserved quantity. The randomness in the examples used in 
this article is introduced through a uniform distribution for 
convenience. However, the conserved area theorem works for 
other methods of randomization such as Gaussian and Poisson 
distributions and is valid even in the absence of ensemble 
averaging. It is plausible that the conserved area theorem can 


be extended to 2D and 3D photonic crystals pl| , which will 
be the subject of future studies. 


VI. Appendix 

Here, we present a proof of the conserved area theorem, i.e. 
fo dfi log (JM^ } |- 2 ) = 27Vlog(r), where is defined 

in Eq.p| Equivalently, we can show dip log = 




0 , where pW is defined as: 


N 


pw=/xm( # )= n 


Wn 


\n= 1 


N 

n 

n= 1 


1 

-ip 


ipe 2llfn 

g 2i<p n 


(9) 

Using Eq. [ 9 ] and defining z = exp(—2 icp), it can be shown 
that 

J dp log(p ^°)=7 + y ~^log(l + x(z,p,{ 0 n}))> 

c (10) 
where 7 — i dp J2 n=1 X is a holomorphic func¬ 
tion of z that vanishes at 2 : = 0; therefore, Cauchy’s 
residue theorem implies that the contour integral vanishes and 
Jo dp log (p 22 ^ = 7 - Therefore, we can show that 

J ‘^r loR ( P'p ' I <Z£log(l 
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( 11 ) 

which completes the proof (c.c. stands for the complex con¬ 
jugate). In the last step, we have used the fact that 7 is a 
purely imaginary number. Note that averaging over {0 n } is 
not needed to prove the theorem; therefore, it is valid for each 
random transmission curve. 
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